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0> : ABSTRACT 

0^ 

'^ , In order to understand the positional uncertainties of arbitrary objects in several 

of the current major databases containing astrometric information, a sample of 
f^ I extragalactic radio sources with precise positions in the International Celestial 

Reference Frame (ICRF) is compared with the available positions of their optical 
CN ■ counterparts. The discrepancies between the radio and various optical positions are 

f~^ ' used to derive empirical uncertainty estimators for the USNO-A2.0, USNO-Al.O, 

C^^ ■ Guide Star Selection System (GSSS) images, and the first and second Digitized Sky 

\^ . Surveys (DSS-I and DSS-II). In addition, an estimate of the uncertainty when the 

^P , USNO-A2.0 catalog is transferred to different image data is provided. These optical 

astrometric frame uncertainties can in some cases be the dominant error term when 

cross-identifying sources at different wavelengths. 
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Subject headings: astrometry 



1. INTRODUCTION 

In the past decade it has become increasingly easier to tie an arbitrary optical field to a 
standard astrometric reference frame with arcsecond precision. With the recent publication of the 
USNO-A2.0 astrometric catalog (Monet et al. 1998) it is now possible to achieve ~0'.'5 accuracy 
at any point on the sky. This greatly enhances the ability to search for optical counterparts of 
sources detected at other wavelengths. This new accuracy is particularly desirable when searching 
for radio sources, which typically have positions more accurate than is achievable for the optical 
frame (see, e.g., Deutsch et al. 1999). In addition, the Chandra X-ray Observatory will for the first 
time provide X-ray positions with arcsecond accuracy, creating a need for better optical precision. 
Indeed, even if the internal precision of a radio or X-ray position is well understood, it is often 
difficult to obtain a quantitative estimate of the additional uncertainty introduced by the transfer 
of these coordinates to a given optical astrometric frame. For many problems, this component can 
be the dominant source of positional error. This work describes an empirical, uniform test of the 
astrometric accuracy of several well known databases containing astrometric information. 



1.1. HST GSC and GSSS Images 

The first astrometric catalog with full and dense sky coverage was the Hubble Space Telescope 
(HST) Guide Star Catalog (GSC), which is described in detail by Lasker et al. (1990). Some 
initial defects in version 1.0 of this catalog have been corrected in the version 1.1 release. The 
GSC contains some 2 x 10^ objects incompletely down to ?ti ~ 15. As it is possible that a deep, 
small field CCD image may contain few, if any, unsaturated GSC stars, it is often necessary to 
refer directly to the Guide Star Selection System (GSSS) images, the plate data used to create the 
GSC itself. In the North, 1982-1985 epoch "Quick V" plates are used, and in the South, deeper 
1974-1984 epoch UK SERC J plates are used. Each plate image is digitized to 1'.'7 pixel"^, and 
the coefficients for a third order polynomial astrometric solution are included in the image header, 
defining the GSC astrometric frame. The astrometric accuracy of the GSSS images is discussed by 
Russell et al. (1990). Although the GSC has been available on CD-ROM for a decade, extractions 
from the GSSS images themselves have only been accessible outside the STScI via the WWW 



( http://archive.stsci.edu/ ) for the last few years. Note that additional astrometric improvements 
are implemented in the GSC 1.2 release (see Roser et al. 1997) but as these improvements are not 
easily incorporated into the GSSS images, this catalog is not considered here. 



1.2. DSS-I 

The next astrometric archive to become widely available was the Digitized Sky Survey 
(DSS-I). This image dataset was distributed on 102 CD-ROMs after undergoing a lossy 10 x 
compression process. In the South, the UK SERC J GSSS images are used, while in the North, a 
digitization of the 1950-1956 epoch Palomar Observatory Sky Survey I (POSS-I) E (red) plates is 
used, again at 1'.'7 pixel"^. Although the astrometric accuracy in the South is identical to that of 
the HST GSC, the astrometric accuracy of the DSS-I in the North is acknowledged to be inferior 
to the GSC. The POSS-I E data were selected for the DSS-I because they are deeper than the 
"Quick V" data, and comparably deep to the southern UK SERC J images. The DSS-I was a 
major milestone as it made a snapshot of the sky down to rra ~ 20 with embedded astrometric 



information available to anyone with access to the CD-ROM set, or shortly thereafter, anyone 
with WWW access to the ST Scl or other sites where the Survey data are available on-line. An 
assessment of the astrometric accuracy of the DSS-I has been previously published by Veron-Cetty 
& Veron (1996), and a more detailed description of the DSS-I itself can also be found therein. 



1.3. DSS-II 

The Palomar Observatory Sky Survey II (POSS-II) and the Second Epoch Survey (SES) 
in the south are currently being processed to form the Second Digitized Sky Survey (DSS-II). 
It contains a 1985-1997 epoch multicolor entire-sky survey using newer, better quality plates 
digitized to slightly more than I'.'O pixel"^. While nearly all of the data have been obtained, only 
about 85% of the sky (red plates only) is available at the ST Scl WWW interface at the time of 
this writing, and it should be further emphasized that these images are being made available with 
prepublication status. The astrometric information assessed here is likely to change when the final 
product is available. Recent status reports of the Digitized POSS-II (DPOSS; does not include 
the SES) have been presented by Djorgovski et al. (1997, 1998). 



1.4. USNO-Al.O 

The USNO-Al.O astrometric catalog (Monet et al. 1996) contains 5 x 10^ objects down to 
m ~ 20, detected on plates digitized by the Naval Observatory's Precision Measuring Machine 
(PMM). For 6 > -33°, POSS-I O and E plates are used, and for 6 < -33°, UK SERC-J and 
European Southern Observatory ESO-R survey plates are used. Only objects which are detected 
in both the red and blue plates (within 1") are included in the catalog, and the catalog positions 
are derived from the blue plates. Both a blue and a red magnitude are also published for each 
object in this catalog, which is available in a 10 CD-ROM set and via the WWW. The absolute 
astrometric frame is based on the HST GSC version 1.1, and is therefore not significantly more 
accurate than the GSC. Whereas the images discussed above come only with a single polynomial 
astrometric solution to describe the entire plate, the USNO catalog positions are based on images 
which are described with a polynomial plus a correction map to compensate for the systematic 
errors inherent at the edges of the Schmidt plates. This allows significantly better precision over 
the entire area of the plates. See also Canzian (1997) for additional information. 



1.5. USNO-A2.0 

The USNO-A2.0 astrometric catalog (Monet et al. 1998) is based on the USNO-Al.O with 
a few important differences. The most important difference is that the absolute astrometric 
calibration is based on the International Celestial Reference Frame (ICRF) as realized by the 
USNO ACT catalog (Urban et al. 1998), which itself is based on the Tycho catalog (ESA 1997), 
generated from the Hipparcos mission. Another important difference is the source of individual 
positions; whereas the Al.O gives blue plate positions, the A2.0 gives the average of the blue 
and red plate positions. POSS-I O and E plates were typically obtained on the same night, but 
SERC-J and ESO-R plates were obtained at different epochs and thus the catalog position will 
place objects with proper motion at the average epoch of the two plates; objects which have moved 



more than 1" between the two epochs will not be included in the catalog at all. The alignment 
with the ICRF makes the USNO-A2.0 the most accurate astrometric catalog of faint stars, i.e. 
the stars likely to be found unsaturated on an arbitrary CCD image. This catalog is available 
on an 11 CD-ROM set and via the WWW from the US Naval Observatory Flagstaff Station 



( http://www.nofs.navy.mil/ ). The image data used, which was digitized at about 0'.'9 pixel ^, 



may also become remotely available in early 2000. 



1.6. GSC-II 

Although not yet available at this time, the Second Guide Star Catalog (GSC-II) will likely 
be available soon, and promises to provide positions comparably accurate to the USNO-A2.0, as 
well as magnitudes, colors and proper motions. This catalog will be based on the DSS-II plate 
material, and the first epoch for the proper motion measurements will be based on the DSS-I 
plates. The GSC-II is merely mentioned here as a future valuable resource for completeness and is 
not discussed further. For additional information on the construction and status of the GSC-II, 
see McLean et al. (1997, 1998). 



2. ANALYSIS 

To assess the astrometric accuracy of these databases, a dataset of objects with extremely 
accurate radio positions is obtained from Ma et al. (1998), who present a careful derivation of 
VLBI positions, most of which have milliarcsecond (mas) accuracy, for 608 extragalactic radio 
sources in the ICRF. These sources are grouped into three classes, of which all are used in this 
work. The handful of sources which have positional uncertainties larger than 30 mas are not used. 
The sources are fairly evenly distributed over the entire sky. The right ascension zero point has 
been aligned to the ICRF which is consistent with the FK5 J2000.0 optical frame. 

For each of these sources, a small 1' x 1' region about each position is extracted from the 
USNO-A2.0 catalog and the DSS-I images. The images are examined individually, and those 
optical counterparts which are not suitable for this study are thrown out. Briefly, 62 sources 
are somewhat subjectively deemed too faint for an accurate centroid based on the DSS-I image 
(not quite the same data as used to create the catalog, however); 64 sources are too crowded by 
neighbors or too extended; 136 sources were too faint and not even in the USNO-A2.0 catalog; 63 
sources were not present in the USNO-A2.0 catalog but yet appeared to have sufficient signal in 
the DSS-I images that they probably ought to have been. This leaves a dataset of 283 sources 
for which the data are acceptable in at least some of the databases discussed here. The datasets 
used for testing of the various databases are sometimes reduced further. For the three imaging 
databases, a minimum of 20,000 ADU ("counts") within a 5" aperture is required; for example 
this reduces the dataset used for the GSSS images by 41 sources as the northern "Quick V" images 
are not as deep as the DSS-I images. Two of the 283 sources are not in the USNO-Al.O catalog 
(in fact those particular fields contain far fewer sources in the Al.O than A2.0 overall). As noted 
before, not all of the DSS-II images are yet available on-line, and thus the accuracy of the DSS-II 
was only assessed with the 235 sources for which a good DSS-II detection was available. 

For the 283 sources, positions from all the available databases, either directly from a catalog 
or centroided on an image, are recorded along with the radio coordinates of Ma et al. (1998) for 



later inspection. The total flux of the sources is also measured using simple aperture photometry 
within a 5" radius. 

In the fields of ~ 100 of the radio sources, a set of isolated, easily measurable reference stars 
with USNO i? < 18 is selected to transfer the USNO-A2.0 astrometric frame to each of the three 
images. The transfer fit usually employs 30 — 80 stars, depending on the field density, and yields an 
uncertainty in the transfer of <0'.'05. The transferred solution is a simple linear solution (sufficient 
for the small 10' x 10' fields used for this purpose) derived with a least-squares fit using procedures 
from the IDL Astronomy User's Library (Landsman 1993) and written in IDL by the author. The 
coordinates of each optical counterpart are then also recorded using the transferred solution. 



3. DISCUSSION 

Using the optical positions of these extragalactic radio sources measured in these various 
ways, the distribution of the deviation from the radio positions is examined. In Table 1 relevant 
derived parameters from these distributions are listed for the different databases. Each astrometric 
database is listed in Column 1. Columns 2-3 list the average Aa and A6, indicating any 
global shift from the ICRF. In columns 4 through 6 the maximum observed Aa, A6, and AR 
(AR = \/Aa^ + A(5^) errors in arcseconds within the best 68 percent of the respective distributions 
are given. The next four sets of 3 columns list similar values within the 90, 95, 99, and 100 
percentiles, where the latter means all of the sources. The number of sources examined for each 
database is provided in the last column. For example, for the USNO-A2.0, 90% of the 283 objects 
examined differed radially from the radio position by less than or equal to 0'.'40. Similarly, 90% of 
the 283 objects (but not necessarily the same 90% as above) differed in 5 by less than or equal to 
0'.'30. These results will be discussed for the various databases in turn. 



3.1. USNO-A2.0 

The USNO-A2.0 is clearly superior to the other frames considered here. The average Aa 
and A6 are quite small but yet deviate from 0.0 in a marginally significant way. The astrometric 
uncertainties are better than the other databases by nearly a factor of 3, with la ~ 0'.'25, 90% 
confidence at 0'.'40, and even 3a ~ O'.'S. In fact, the largest observed discrepancy in this sample of 
283 sources is 0'.'63! This is quite remarkable given the possibility of true misalignment between a 
radio source and a slightly extended optical galaxy, not readily discernible as extended on the low 
spatial resolution DSS-I images. A brief examination of this most discrepant source, 0458-020, 
yields no anomalies in the images and is cataloged as a high-z QSO, and thus not likely to suffer 
from intrinsic radio/optical misalignment. Centroiding errors probably also contribute to these 
overall uncertainties in addition to local frame offset errors, but it is the combined effect that it 
most relevant in most cases. 

In Fig. 1 the distribution of discrepancies is plotted against a variety of other parameters. 
The first panel depicts the discrepancies in Aa and AS in arcseconds, in the sense that the radio 
coordinate is subtracted from the optical coordinate. The distribution is quite tight and slightly 
offset from the center. A fiducial 0'.'5 radius circle is overplotted. The second panel shows the 
combined (radial) discrepancy (AR) vs. Declination; no Declination region is apparently worse 
than others, even near the poles. The third panel shows AR vs. the distance of the source from 



the plate center on the DSS-I images, which is mostly the same plate material as used to create 
the USNO-A2.0. The units are in 1'.'7 pixels as encoded in the image headers. The horizontal 
solid line is simply the O'.'S fiducial. The vertical dashed line indicates the approximate useful half 
width of the plates. Clearly the astrometric solution is essentially as good near the edges of the 
plate as near the center. 

The fourth panel shows AR vs. the USNO R magnitude. It does appear that the faintest 
stars post slightly worse values. If all stars i? > 18 are excluded, the resulting errors are ~ 10% 
better for the new sample of 225 objects (see Table 1). In the fifth panel AR is plotted vs. the 
USNO {B — R) color. Even the bluest and reddest objects show no obviously greater errors than 
the others. Finally in the last panel AR is plotted vs. log^Q ADU within a 5" aperture around the 
source in the DSS-I image; this is not directly indicative of the brightness in the images used to 
create the A2.0, but a reasonably good estimate of it. A similar figure is generated for each of the 
databases examined. They are not all reproduced here, but are available from the author. 

Figure 2 displays the distributions of Aa and A5 in arcseconds for the USNO-A2.0 283 object 
sample. The slight average offsets evident in Table 1 can also be readily seen here and appear 
to be significant. The reason for this systematic offset remains unclear. Certainly for program 
objects of color significantly different from the typical USNO-A2.0 star, atmospheric differential 
refraction must play some role. 



3.2. USNO-Al.O 

The USNO-Al.O catalog is worse than the A2.0 by nearly a factor of 3, and shows similar 
results as the GSSS images, not surprising as the Al.O is based on the absolute astrometric 
calibration of the GSSS/GSC. The results are probably slightly better due to improved plate 
distortion corrections. A set of plots like Fig. 1 for this database turns up nothing unusual. 



3.3. GSSS Images 

The GSSS images provide the next best astrometric database. It is still in use for HST 
observation planning. Russell et al. (1990) assessed the astrometric accuracy of the GSC and GSSS 
images. Using a small sample of 48 extragalactic radio sources, they derived aa = 0'.'63, as = 0'.'58, 
although the values were better in the north. However, with their much larger "CAMC" sample, 
they find for the whole plate aa = 0'.'58, as = 0'.'53, and slightly better, Uq, = 0'.'55, as = 0'.'48, for 
the inner 50% of plates. This compares favorably to the aa = 0'.'47, as = 0'.'47, ur = 0'.'71 values 
derived here, and lends considerable confidence that the reduction described herein is accurate. 



3.4. DSS-I 

The DSS-I begins with a reasonably good la astrometric uncertainty but degrades at > 99% 
confidence to the worst statistics of the group. The worst offset (for the object 0202-1-319) is 4'.'6. 
The accuracy in the South is identical to that of the GSSS images, as they are the same data. In 
the North, the POSS E images have inferior astrometric fit quality, as is noted in the documents 
which come with the CD-ROMs, although it is perhaps not widely appreciated how far off the 



astrometry can be. A figure like Fig. 1 for this dataset indicates that the worst discrepancies he 
near the outer edges of the plates. If all sources greater than 4500 pixels (2? 125) from the plate 
centers are excluded, the sample of 144 remaining objects shows a somewhat better distribution as 
seen in Table 1. Clearly, positions derived from the DSS-I should be treated with caution as there 
is a potential to be several arcseconds in error in the North. A far more detailed discussion of 
the DSS-I astrometric accuracy is presented by Veron-Cetty & Veron (1996). They derive similar 
uncertainties as in Table 1 and derive some corrections. 



3.5. DSS-II 

Immediately noticeable in the DSS-II row in Table 1 are the large <Aa>, <A5>. In the plot 
like Fig. la for the DSS-II it is immediately evident that most sources exhibit a (+1", —1") offset. 
As the pixelsize is ~ 1" it is likely that this is due to a one pixel offset somewhere in the pipeline 
which generates the astrometric solution for these images. It is, however, also readily evident that 
a smaller subset of discrepancies do lie about (0,0) and thus some images may not suffer from this 
defect. For this reason, it may not be prudent to simply apply this offset to the images. A second 
set of statistics is generated for the DSS-II with all source (a, 6) positions adjusted by (—1", +1")- 
The result is better, but still not as good as the GSSS image or USNO catalogs. It should be 
stressed that these results are derived from the prepublication data available at this time, and the 
astrometry in the final product may be significantly improved. 



3.6. Transferring the USNO-A2.0 Frame to Images 

The uncertainties involved in obtaining an ICRF position for an extragalactic point source 
from the USNO-A2.0 catalog is now fairly well established. The positions of objects in the Galaxy 
with known proper motion can also readily determined, although determining the exact epoch 
of the A2.0 positions is nontrivial (the epoch for each field is available in a separate catalog at 
the NOFS WWW site). For objects with unknown proper motion, the uncertainties must be 
estimated on a case by case basis, depending on the data available. 

The remaining case is for objects simply not listed in the USNO-A2.0 catalog, for whatever 
reason. Since large telescopes with modern detectors can go much fainter than the A2.0 is a 
relatively short time, this is actually a frequent issue. Obtaining coordinates for such objects is 
easily done by fitting an astrometric solution to the new image using an input list of USNO-A2.0 
(a, 5) and image (X, Y) positions for a set of stars common to both the catalog and the image. 
Given a few tens of reference stars, the final uncertainty in the transferred solution due to 
centroiding errors can easily be less than 0'.'05. However, if the epoch of the USNO-A2.0 and 
image data differ significantly, the small proper motions of the reference stars will dominate 
the uncertainty. Stars with large proper motions can, of course, be easily rejected from the fit. 
If the directions of motions of the reference stars are mostly random, the residuals are merely 
inflated. The worst scenario is where the tiny proper motions of many reference stars yield a 
small systematic shift due to Galactic rotation. This can be beneficial for stars moving with the 
reference stars, but is detrimental to determining positions of extragalactic objects, or objects 
otherwise not associated with the bulk motion. In an effort to quantify the uncertainties in 
determining positions of extragalactic objects not available in the USNO-A2.0, an A2.0 solution 



is transferred to several of the database images used here and the resulting position of the radio 
source counterpart measured as described in §2. The results for the DSS-II images may be most 
relevant as they have fairly recent epochs and are better sampled and deeper than the other 
images used here. 

The results of this effort are tabulated in the last three rows of Table 1, where column 1 lists 
the source database which is aligned to the USNO-A2.0. The results yield a ~ 25% degradation 
over the A2.0 itself, which is quite respectable given that the difference in epoch between A2.0 
images and DSS-II images can be up to 40 years. For many applications these larger values will be 
a better estimate of the true uncertainty of measured positions based on, but not directly from, 
the USNO-A2.0. 

In closing the reader is reminded that many precise positions are still found in B1950, and 
special care must be taken when converting them to J2000 if subarcsecond accuracy is important. 
Some generic precession utilities do not perform the necessary FK4 to FK5 correction, and the 
resulting improperly converted positions may be in error by ~ 1" . 
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Table 1. Empirical Uncertainty Estimates in Arcseconds Relative to the ICRF 





<Aa> 


<AS> 




68% 




Ac 


90% 
AS 


AR, 




95% 




Aq, 


99% 
A,5 


AR, 




100% 






Database 


Ac 


AS 
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0.04 


0.04 


0.17 
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0.26 


0.29 
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0.40 
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0.34 


0.45 


0.45 


0.43 


0.52 
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0.59 


0.63 
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0.03 


0.16 


0.17 
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0.26 
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0.34 
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0.43 
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0.45 


0.45 


0.53 


225 
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-0.11 


-0.12 


0.47 


0.41 


0.65 


0.78 
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1.04 


0.94 
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1.41 


1.70 
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-0.12 


-0.07 


0.47 
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1.31 
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4.59 
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0.81 
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1.37 
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1.68 
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0.74 
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Fig. 1. — The distribution of discrepancies between USNO-A2.0 and VLBI radio positions (Ma et al. 
1998) for the sample of 283 extragalactic radio sources. A fiducial O'.'S radius circle or horizontal 
line is overplotted in each panel, a) Act vs. A6 in arcseconds; b) AR (radial discrepancy) in 
arcseconds vs. Declination; c) AR in arcseconds vs. distance from the plate center (1'.'7 pixels) 
where the dashed vertical line denotes the plate half width; d) AR in arcseconds vs. USNO-A2.0 R 
magnitude; e) AR in arcseconds vs. USNO-A2.0 {B — R) color; f) AR in arcseconds vs. log^^g ^DU 
within a 5" radius aperture on the DSS-I images. There do not appear to be any large systematic 
trends in any of these representations. 
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Fig. 2. — Histograms of Aa and AS discrepancies between USNO-A2.0 and radio positions for a 
sample of 283 extragalactic radio sources. There do appear to be slight zeropoint offsets in the two 
distributions, although their origin and significance remains unclear. Certainly color terms must 
play some role. 



